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Abstract 

We consider the problem of calculating the Green's functions associated to a massive scalar field 
with modified dispersion relations. We analyze the case when dispersion is modified by higher 
derivative spatial operators acting on the field orthogonally to a preferred direction, determined by 
a unit time-like vector field. By assuming that the integral curves of the vector field are geodesies, 
we expand the modified Klein-Gordon equation in Fermi normal coordinates. By means of a 
Fourier transform, we find a series representation in momentum-space of the Green's functions. 
The coefficients of the series are geometrical terms containing combinations of the Ricci tensor 
and the vector field, as expected from previous calculations with different methods and for specific 
backgrounds. 

PACS numbers: 04.60.m 
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I. INTRODUCTION 



In recent times, there has been a considerable amount of investigations on the trans- 
Planckian problem. As a satisfactory theory of quantum gravity is not available yet, some 
researchers have focussed on the possibility that the unknown physics beyond the Plank scale 
appears in the form of a modification of the dispersion relations above a certain energy-scale. 
The consequences of this assumption have been widely examined in relation to the inflation- 
ary cosmology , the Hawking radiation |3|, and the Unruh effect |4|]. In some cases, 
modified dispersion relations (MDR) locally break the Lorentz-invariance (for a review, see 
[5I). However, one can keep the general covariance of the action, by assuming the existence 
of a preferred frame, encoded by a unit time-like vector field, which breaks dynamically 
the Lorentz invariance In this setup, often named "aether theory", the vector field is 
independent of the metric, in the sense that it represents a new degree of freedom in the 
Einstein-Hilbert action. In the most general case, the action contains several new terms as- 
sociated to the vector field, its covariant derivatives and contractions with the Ricci tensor. 
All these terms carry coefficients which are constrained by observations, and recent results 
show that the aether theory is not incompatible with available data (see, e.g., 

The vector field can be coupled to covariant derivatives of quantum fields, thus generating 
modified dispersion relations. In flat space, the propagator associated to these fields contains 
higher powers of the spatial momenta, which displays explicitly the breaking of the Lorentz 
group down to a spatial rotation subgroup. Thus, also in curved space one expects that the 
Hadamard form of the propagator is modified. As a consequence, the fundamental structure 
of quantum observables is certainly different with respect to the Lorentz-invariant case. As, 
in the semi-classical theory, the expectation value of the renormalized stress tensor is the 
source term of the Einstein equations (see for instance Q]), the calculation of the modified 
Green's functions is of fundamental importance in order to establish the quantum back- 
reaction, and, hence, the effects of trans-Planckian modes. Some investigation on this issue 
has recently begun, especially in the context of cosmology. For example, the renormalized 
stress tensor has been calculated, via adiabatic regularization, in the cases of FLRW and 
Biaachi I ^logical ba„ s fl Q. How™, it wouid be desirable to bave a g ene r ai 
expression, which can be applied to any background. 

The aim of this paper is to find a general expression for the Green's functions with 
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modified dispersion relations. In a previous paper we found a representation of these 
functions in momentum-space, in the case when the background is ultra-static. Here, we 
generalize this result for any background and for any analytic (in momentum-space) disper- 
sion relation. When the latter is Lorentz-invariant, there are various equivalent techniques 
to calculate the Green's functions, such as the deWitt-Schwinger expansion in coordinate 
space 12j, and the Bunch and Parker expansion in momentum-space 13J]. When MDR are 



present, however, the first technique is difficult to apply as it heavily relies upon Lorentz- 
invariance, while the second does not need this symmetry. In fact, Lorentz-invariance in the 
Bunch and Parker expansion allows to reduce the number of terms in the expansion, but it 
is not essential, as already shown for the ultra-static case in JjJ. For general backgrounds, 
however, the Bunch and Parker method leads to algebraically complicated results because 
of the Riemann normal coordinates (RNC) used in the expansion of the operators. This 
coordinate system is very useful when one needs to expand the metric around one fixed 
point. But when, besides the metric, there exists also an independent vector field, together 
with differential operators associated to it, the expansion of the latter becomes very compli- 
cated in RNC. As the vector field define a congruence of time-like curves, it would be more 
convenient to find a coordinate system adapted to these. 

Such a system exists and it was discovered by Fermi Jjj], and developed later by several 
authors (see, e.g. Q)- It generalizes the RNC, in the sense that it allows to define a 
local orthonormal coordinate system non only about one fixed point, but about an entire 
geodesic curve \y\ . Thus, the use of the so-called Fermi normal coordinate (FNC) allows to 
simplify the problem of expanding the operators associated to the vector field, and provides 
for a well-defined orthonormal frame. With this frame, as we will show, it is possible to 
construct a momentum-space representation of the Green's function by Fourier transforming 
the Klein-Gordon operator. 

The plan of the paper is the following. In Section II, we introduce the modified Klein- 
Gordon operator in the context of the aether theory. In Section III, we calculate the 
momentum-space representation of the Green's functions. Most of the technicalities are 
confined in the appendices: in the first, we briefly sketch the FNC construction, while, in 
the second, we display the expansion in FNC of the differential operators contained in the 
Klein-Gordon equation. We conclude, in Section IV, with a discussion of the results and 
few remarks. In this paper, we set h = c = 1, and the signature is (—,+,+,...). 



3 



II. AETHER THEORY AND MDR 



In the aether theory, one considers the Einstein-Hilbert action implemented by terms de- 
pending on a unit time-like vector field u^. The general form in N dimensions is given 
by 

5 = T7T7^ / d N x^9(R-2A + C u ) , (1) 

l07TG J 

where R, A, and G are the Ricci scalar, the cosmological constant and the Newton constant 
respectively. The term C u is the Lagrangian associated to the vector field, and its general 
form reads 6] 

C a = -A(u% + 1) - hF^F^ - hiyXf - hR^u v - b A u a u p V a u„V pu" . (2) 

In this expression, F^ = 2V[ M Wj,] , and the constant A acts as a Lagrange multiplier, which 
ensures that the vector field is unit and time-like. Finally, is the Ricci tensor, and 
the coefficients b\ . . . 64 are arbitrary constant. The form of this Lagrangian guarantees 
that the theory is generally covariant, even though the vector field dynamically breaks local 
Lorentz invariance, as it determines a preferred direction. Some of the terms in Eq. (J2J) 
vanish whenever the integral curves of the vector field are geodesies. In such a case, in fact, 
the vector field is everywhere hypersurface orthogonal and F^ = 0. Also, the acceleration 
a /j. = u «y aU n = q ? anc [ the term proportional to 64 vanishes. In the following, we assume 
that the integral curves are geodesies, thus the observer comoving with w M is free-falling. In 
any case, to analyze the equations of motion, it is always convenient to choose a background 
of the form 

g^dx^dx" = -(u a dx a ) 2 + q^dx^dx" , (3) 

where = g^ v + u^u v is the induced metric on the hypersurfaces orthogonal to j^, E|. 
We now consider a free, massive, minimally coupled scalar field with a Lagrangian given by 



r _ 1 



(4) 



where the operator V is defined such that 

V 2 <P = q\V a (q^V^) , (5) 



and the coefficients A pq are arbitrary {^.f]]. This Lagrangian explicitly breaks the local 
Lorentz group to a rotational subgroup 6(. The Klein-Gordon operator associated to the 
field can be written in the form 



(□ - m 2 )(j) + F [V 2 ] = 



(6) 



where F is some functional of the operator V 2 , such that it can formally be expanded as 

oo 

F[V 2 ]<f ) = Y,^nV 2 y, (7) 



n=2 



with arbitrary coefficients a^n- This operator is responsible for the modification of the 
dispersion relation as, in momentum-space, it introduces higher powers of p 2 , the square of 
the spatial momentum. For our purposes, it is more convenient to write D0 as 



□0 = V A <p - u a u p V a V P (t) - K u a d a ( 1 



(8) 



where K = q^ u ^7 ^u v is the trace of the extrinsic curvature of the hypersurface orthogonal 
to m m . In this way, we can write the Klein-Gordon equation as 



V 2 <p - u a u p V a V P (j) - K u a d a cf> - m 2 4> + F [V 2 ] = 



(9) 



If F can be expanded as in Eq. (CO), we can include the first term of the above equation into 
the last one, and finally write the iV-dimensional Klein-Gordon equation as 



«2n£> 2 > - u a vPV a Vp<t> - Ku a d a (p - m 2 (f) = 



(10) 



n=l 



with the requirement that «2 = 1, while the other a's are arbitrary. 



III. MOMENTUM-SPACE REPRESENTATION OF THE GREEN'S FUNCTIONS 



In this section, we find an expansion up to second order of the metric in FNC, of the Green's 
functions G(x, x') associated to the scalar field <fi of Eq. (flUI) . These must satisfy the equation 



a2n£> 2n - uVV a V p - K u a d a - 



m 



71=1 



G(x,x') 



1 

V9 



8(x — x) 



(11) 



where x and x' are the coordinates of two nearby points, and the differential operators act, 
by convention, on x. In order to write this equation in terms of symmetric biscalars, we 



need to define Q, Q 



G(x,x') = g{x)- x l A G{x,x')g{x')- l,A 



fj 



A/2 (x)5(x-x') = (T 1/4 (x)<5(x - x')g- l/ \x') . 



(12) 



If we choose the FNC as local coordinate system, we can identify P as the point with 
coordinate x' — (0, 0, . . . , 0), and Q as the point with coordinate x = (r, x a ), see Fig. ([T]) in 
Appendix I. Hence, the Green's functions G must satisfy the equation 
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1/4 



a 2n V 2n - u A u B V A V B - Ku A d A 



n=l 



(g- 1/4 G) - m 2 G = -5(x a )5(r) 



(13) 



The convention is that the indices (A, B, . . .) label space-time coordinates while (a, &,...) 
spatial coordinates only, see Appendix I. We now wish to expand this equation in FNC up 
to the second order. By using Eqs. (lAlOh and (1B4[) . we find first the expansion in FNC of 



g^V^g-^G) = 5 ab d a d b G + hiG + Q A b x b d A G 



where we defined 



-U b , V b — -K 06 



H = 6 ab H ab . 



By recursion, it is easy to show that 



g^V^ig'^G) = [d 2n + -Hd 2( - n ~V + nQ A b x b d A d 2 ^ + n(n - l)Q\ d A d b d 2 ^ 



(14) 

(15) 

G, 
(16) 



where we use the convention d 2 G = 5 ab d a dbG, and d k G = G for k = . With Eqs. ( 1A10I) . 
( 1B3I) . and (1A14|) we find the expansion of the remaining terms, and the expansion of Eq. 
( TT3|) in FNC and up to two derivatives of the metric finally reads 



«2n \d 2n G + -Hd 2(n - l) G + nQ A b x b d A d^ n -^G + n{n - l)Q A b d A d b d 2{n ~ 2) G 



n=l 



+ 



- (1 + QabX a x b ) d d Q G - Q a b x b d a G + 3Q° c x c d G - m 2 G = -5(r)8(x a 



(17) 



The complexity of this equation can be greatly reduced by performing the iV-dimensional 
Fourier transform 

d N k 



G{x,x') 



(2k) 



JgABk A x t 



g(pIp 2 " 



N 



(18) 



where k A = (po, p) is the A^-momentum, and x B = (r, x) are the FNC. The arguments of G 
stress the rotational invariance assumed before, together with time-reversal symmetry. As 
discussed in Appendix EJ the geometrical coefficients are considered r-independent at any 
given order, thus the Fourier transform can be calculated to give, after some algebra, 



(S + p 2 -m 2 ) G + 



—-HG + Q a a G + Q \pa9 G 



DS + Q \paP GD 2 S + 



QabP 2 d a d b G + Q\G + Q\p a d b G - 3Q" bPo d b G = -1 . 



(19) 



In this equation, we have defined 



s = j> 2n (-i)y n , d a c 



n=l 



— DS= — 

dp a ' dp 2 



(20) 



To solve this equation, we proceed by iteration, similarly as in [13(. We first write G 
Go + G 2 , where 

1 



Go 



(21) 



m 2 — S — Pq ' 

is the propagator in flat space. Then, we replace G with Go + (?2 i n Eq. (fl9l) and we keep 
only second order terms. Finally, by using the definition ( |2T1) to write S in function of Go, 
we find 



G 2 = - l -HDGo + Q\poP b D{D-?,Do)Go + Q a a (D + D )G - p 2 DD G 



+ 



+ QabPy 



D(D + Do)Go-4p 2 GoD 2 Go 



(22) 



where DqG = dG/dp^. This is our main result, and it displays the corrections to the flat 
propagator, Go, in momentum-space, in terms of geometrical coefficients depending on the 
background, and operators acting on G - These coefficients contain up to two derivatives of 
the metric. Note that, in this form, the equation is independent of the dispersion relation 
adopted. When the dispersion is linear, S = —p 2 and DG = —D Go = — Gq. Thus, the 
above equation reduces to 



Gf 1} = -HGl + 2Q a aP 2 G d + 8Q 0b pVGl - 8Q abP a pY Gt ■ 



(23) 
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IV. DISCUSSION 



We now discuss our results, beginning from Eq. (1231) . This equation should be contrasted 
with the well-know Bunch and Parker expansion [l3|, which, to second order, just reads 



Gf l) = l -RGl, G = (p 2 -pl + mY 1 , (24) 

where R is the Ricci scalar. Our expression seems much more complicate, and the reason 
is due to the choice of FNC instead of RNC Thus, to transform Eq. into the expres- 
sion above, we should go back to Eq. ffTUl) and expand in RNC, together with imposing 
Lorentz-invariance instead of rotational invariance, as we did with Eqs. (1B5|) . In opposition, 
by using FNC, even with Lorentz-invariant dispersion, the expansion is obtained around the 
Minkowski metric and a preferred direction with = <5q at the lowest order. However, 
by integrating both expansions in po, one can see that, in the coincidence limit, the diver- 
gent contributions in Eq^ fl23|) add up and coincide with the divergent part of the Bunch 
and Parker expansion [21[. Therefore, the fundamental Hadamard structure of the two 
propagators is the same 

For Lorentz-breaking dispersions, we have found Eq. (122|) . The geometrical coefficients 
depends on various components of the Riemann tensor. Some of these terms become trivial 
if one assumes that the background has some symmetry. By using the usual transformation 
law of tensors under coordinate change, R a bcd = Rap-yS^a e 7 c e S d, it is immediate to show 
that, in the case of ultra-static metrics, only the first term survives, and we recover the 
results of Ref. [11]. 



In the case of a Bianchi I cosmology, the term Q fe vanishes, while the remaining ones, 
when transformed into comoving coordinates, give rise to terms proportional to R and 
R fiu u fl u u . These terms were also found via adiabatic regularization in Ref. [lo| . however a 



detailed comparison involves the calculation of the inverse Fourier 



transform of Eq. ff22l) . 



• 111 



even in the simple 



This appears quite difficult to achieve. In fact, as we showed in Ref. 
case of a quartic dispersion and ultra-static metric, a coordinate space representation of 
the Green's function is hard to find by inverse Fourier transform of the momentum-space 
representation. 

The future plan is to use our results to regularize the UV divergences of the two-point 
function (and of the stress-energy tensor) with MDR, by operating solely in the momentum- 
space. This procedure is essential to renormalize the stress tensor in the case of static 
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space-times, such as black holes. For these backgrounds, the Klein-Gordon operator con- 
tains higher-order spatial derivatives and WKB methods do not apply, in contrast to the 
cosmological case. Work in this direction is in progress. 

Finally, we remark that our calculations can be extended to more general cases. For 
instance, the coefficients of Eq. (I22p contain up to two derivatives of the metric. Higher 
order expansions are of course possible, although they are likely to be algebraically rather 
complicate. Another extension can be obtained by considering a non-geodesics vector field. 
In this case, one can make use of the generalized Fermi coordinate together with the Fermi- 
Walker transport, which replaces the usual parallel transport and keeps in account the 
acceleration of the local frame Q] • 
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APPENDIX A: FERMI NORMAL COORDINATES 



In this appendix, we briefly review the construction of the Fermi Normal Coordinates (FNC). 



More details can be found in Refs. 



15]- 19] . Let 7 denote the integral curve of the vector 



which corresponds to the velocity of the observer. We assume that the latter is free-falling, 
so that 7 is geodesies. The curve is parameterised by r, the observer's proper time. Let P 
be a reference point on 7, where we set Tp = 0, and Q a point off the curve. The latter is 
chosen to be close enough to 7, so that there exist a unique geodesies (3 orthogonal to 7 an 
passing through Q. We denote by R the point of intersection between 7 and (3. This point 
is located at r = tr. Finally, the vector tangent to (3 at R is called see Fig. (Tj[|). At the 
point P we construct an orthonormal tetrad {e^} such that = u^. We use the convention 
such that frame indices are labeled by A = 0,1,2,... or a = 1,2,..., while general space- 
time indices are labeled by Greek letters. Since 7 is a geodesies, we can construct the same 
tetrad anywhere on the curve by parallel transport. Thus, at any point on the curve, the 
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FIG. 1: Construction of Fermi coordinates. 

components of the tetrad satisfy 

g^u v = -1 , g^ e£ u v = , g^ e£ e v b = S ab . (Al) 

The dual tetrad components are defined by 

el = -u,, e«=5 ab g^. (A2) 

With these, we have the following completeness relations 

<T = + g,u = -ele° l/ + 5 ab e;e b u . (A3) 

The FNC coordinates of the point Q are defined as 

t q = t r , x a Q = -e a a (R)a a (Q,R) , (A4) 

where a a (Q,R) = — £ a is the gradient of the Synge world function relative to (3. By using 
the completeness relations, one can show that 

8 ab x a x b = 2a(Q,R) = s 2 , (A5) 

where s is the spatial distance between the points Q and R along the geodesic (3. By 
expanding the function a in s, one finds the local metric around 7, which has the form 

ds 2 = (r}AB + h AB ) dx A dx B , (A6) 
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where t]ab is the Minkowski metric and 

h 00 = -R 0c0d x c x d + 0(s 3 ) , 



(A7) 



h 0a = ~R cadX C X d + O(s 3 ) , (A8) 

h ab = - l -R acbd x c x d + 0{s 3 ) . (A9) 

These expressions represent the expansion to second order of the metric tensor around the 
curve 7 [22]. Note that indices are raised or lowered with the frame Minkowski metric t\ab- 
The components of the Riemann tensor appearing in these expressions depend in general on 
tr, as they are evaluated at the point R. If the point Q is on the curve 7, the above metric 
reduces to the Minkowski one. The determinant of the metric can be computed by the 
identity det(J + A) ~ 1 + Tr(A), as the tensor Uab can be considered as a small correction 
to the Minkowski metric. Thus, the modulus of the determinant reads 

det \g AB \ = g = 1 - (r c0(1 + ^R a ca ^j x c x d + 0(s 3 ) = 1 - H cd x c x d . (A10) 

In the theory discussed in this paper, besides the metric tensor, there is also the vector 
field By construction, this vector is unit time-like and it is the component of the tetrad 
tangent 7. However, when the tetrad is parallel-transported along f3 to a point off the curve, 
the vector w M shows non-trivial components along the space-like directions. Suppose that 
these can be expanded in a Taylor series around the point R as 

u A = u A \ R + d b u A \ R x b + ^d c d b u A \ R x b x c + ... , (All) 

with u A \ R = 5 A . Then, parallel transport along ft implies £, b V b u A = for all £ 6 . As the 
connection coefficients vanish along 7, in opposition to their derivatives, it follows that 

8 b u A \ R = -T A bD u D \ R = , (A12) 



d c d b u A \ R = -d c Y A bD u u \ R = -R A bc0 , (A13) 

as it results from the expansion of the connection coefficients, which can be found directly 
with the help of Eqs. flA71)- (TA~9]) . and read 

A 2 A AAA 

r afe = -3 R ( ab ) c x c , r os = T B0 = R Bc0 x c . (A 14) 
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Thus, the second term in the expansion (1A11I) vanishes, and we find 



A -5 A -^R A ij0 x i x j + O(s 3 ) . (A15) 

With these components we can finally calculate the expansion, to second order, of the 
induced metric, which reads 

q AB = ^AB + S A S B _ h AB _ Sq (A R B)^ j j + Q ^ (Alg) 

APPENDIX B: DIFFERENTIAL OPERATORS IN FNC 

In this appendix, we find the expansion in FNC of the differential operator V 2 , which appears 
in Eq. ffT0]) . The quickest way is to use the definition {2] 

V 2 f = q AB V A V B f + Ku A V A f , (Bl) 

where K is the trace of the intrinsic curvature of the hypersurface orthogonal to u M . As 
Kab = q c Aq D b^c'Ud, with the results of the previous appendix, we find that the only 
non- vanishing terms are 

K ab = ^R obac x c . (B2) 

We notice immediately that the extrinsic curvature is not symmetric. This is a result of the 
fact that the vector u A is hypersurface orthogonal only at the point R, see Fig. ([I]). However, 
when we parallel-transport the Fermi tetrad to the point Q along /?, orthogonality is lost, 
and this gives rise to a non-zero vorticity at Q, which, in turn, renders Kab n °t symmetric 
[20! ] . For our calculations however, we just need the trace of Kab-, which is given by 

K = q AB V A u B = -\R\ac x c = —Roc x c . (B3) 

By putting together these results, we finally find 

V 2 f = 5 ab d a d b f - \R ab ac x c d b f - l -R\ x a d f . (B4) 

In order to arrive to this expression, we made use of the symmetries of the Riemann tensor 
and the rotational invariance of the propagator. In fact, if / is rotationally invariant (i.e. it 
depends on powers of x 2 ), we have 

h°%d b f = ~R° c b d x c x d d d b f = . (B5) 
12 



and 

h ab d a d b f + ^R ad ac x c d d f = 0. (B6) 

We conclude this appendix with an important remark. We mentioned before that the co- 
efficients of the expansion of g AB depends on r. Therefore, when we compute terms like 
g l 1 4 T> 2 (g~ l 1 4 f) , as in Sec. Ill, we need to know, in principle, the explicit form of the r- 
dependence of the Riemann tensor. In fact, there are r derivatives acting on g, defined by 
Eq. (1A10I) . Since we are interested in the short-distance behaviour of the Green's functions, 
we can take Q close to P, which in turn leads to r R close to t p , see Fig. (CQ). Hence, we 
can expand in Taylor series the coefficients Ro c od, Rocab, and R ac db with respect to r. How- 
ever, the results would be coefficients containing more than two derivatives of the metric. 
Thus, in general, if we fix the order of the expansion to n, say, then any r-derivative of the 
coefficients yields terms of order n + 1 or higher, which can be neglected. This means that 
the coefficients appearing in Eqs. (lA7j) - (lA9j) can be consistently considered as constant in 
r when hit by a r-derivative. A further consistency proof follows from the fact that the 
operator □ can be defined as 

□/ = g~ 1/2 d A (g 1/2 g AB d B f) , (B7) 

where r-derivatives of g explicitly appear, or as 

□/ = g AB v A v B f . (B8) 

The two expressions coincide if, and only if, the coefficients are treated as time- independent. 
For completeness, we report the expression of □/, which reads 

□/ = ( V AB - h AB ) d A d B f - ~ R Ab db x d d A f - R ob od x d d b f . (B9) 
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